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A plasma was produced by the interaction of a 300 MW, 25 nsec, Nd 
laser pulse with an aluminum target. The resulting plasma expanded into 
an ambient background of 2.5 X 10”^ torn and was analyzed using floating 
double probes, magnetic probes and a quartz pressure probe. An early 
disturbance was noted before arrival of the main plasma. 

Further experiments separate this early signal into photoelectric 
response and two fast plasma pulses traveling with constant speeds of 
1.1 X 10^ cm/sec and 5-9 X 10 ^ cm/sec. 

Mapping of the plasma density indicates that the early plasma pulse 
is not syrrmetric with respect to the target normal but expands along a 
line defining the reflected laser pulse. This same mapping indicates 
that the main plasma expands anisotropically for the first 60 nsec 
resulting in an early time asymmetry with respect to the target normal. 
At tines greater than 120 nsec, the asymmetry of the main plasma is no 
longer evident and its density distribution is symmetric with respect 
to the target normal. 
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I. INTRODUCTION 



Early investigations concerning laser induced electron emission by 
William I. Linlor [Ref. 36] in 1962 were followed by a multitude of 
papers dealing with laser produced plasmas. Lichtrnan and Ready [Ref. 

3*0 were among the first to detect an early signal preceding the main 
plasma body. They ascribed the signal to reverse photoelectric effect 
[Ref. 3^] from the electrostatic probes being used. This hypothesis 
received immediate criticism from Honig [Ref. 23] who ascribed the 
early signal to a fast thermal effect. Since then, several new 
hypotheses including ionizing potential waves [Refs. 2 and 31] and 
detonation or blast waves [Refs. 9 and 22] have been forwarded to 
explain this early disturbance. However, a review of the literature 
does not reveal a conclusive report on the behavior or content of this 
signal. 

The separate discovery of spontaneously generated magnetic fields 
by Korobkin and Serov [Ref. 29] and by Stamper [Ref. 43] further 
stimulated interest in this area. Computer simulations have indicated 
the presence of megagauss magnetic fields early in the plasma expansion. 
Fields of this magnitude will influence plasma containment schemes 
generated in the controlled fusion effort. 

In response to an obvious need for complete understanding of the 
laser produced plasma and its induced magnetic fields McKee [Ref. 38] 
and Bird [Ref. 8] undertook extensive surveys of the self-induced 
magnetic fields and attempted to relate them to the expanding plasma. 
During Bird's investigation he noted the presence of an early magnetic 
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signal. This revelation led to renewed interest in the early electro- 
static probe sigials. There are four possible sources for this early 
sigjial . They are : 

1. laser induced electron emission. 

2. East plasma blowoff . 

3. Electromagnetic disturbances. 

i +. Ionization of the background gas by plasma produced photons or 
ionizing potential waves. 

This investigation has been conducted in an effort to determine which 
of these causes is responsible for the observed signal and to provide 
a more detailed understanding of the behavior of the phenomenon producing 
the signal. 

The thesis is divided into six more sections and two appendixes. 
Section II is a description of the general experimental arrangement . 
Section III describes the experiments conducted to solve the problem 
and discusses the results. Section TV attempts to relate two of the 
observations made in Section III to theory. Section V presents a 
summary with conclusions. Section VI provides recommendations for 
further investigation. Appendixes A and B describe the construction 
and theory behind two of the diagnostic tools used during this research. 
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II. EXPERIMENTAL ARRANGEMENT 



A. LASER 



The laser used was a Korad K-1500, Q-switched neodymlum-in-glass 
system. Ihe system consists of an oscillator which generates the 
initial 25 nsec (full width at half power) pulse through Q-switching by 
Pockels Cell. The oscillator pulse is expanded and passes through an 
amplifier which increases the energy of the pulse. The arrangement is 
shown in Figure 1. A detailed description is provided by Davis [Ref. 7]* 



CW Alignment Laser 

i 

Pockels Cell 



Polarizer — 



<? 

0 



K-l Nd Oscillator 





Figure 1. Block Diagram of the General Experimental Layout 
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The output of the laser was controlled by varying the voltage 
applied to the K-l and K-2 flats hi amps. The range of energies 
obtainable is 2.5 to 12 Joules. The duration of the pulse is 25 nsec 
giving a power range of 100 to 480 MW. The laser focal spot size is 
approximately 0.15 cm in diameter after focusing. mu -~ — • — - ^ 



B. LASER MONITORING TECHNIQUES 

Approximately one percent of the laser pulse energy was reflected 
by a beam splitter onto an MgO block. The laser output was monitored 
from this block with a Korad K-Dl photodiode after passing through a 
0.1 percent neutral density filter. This photodiode provides one signal 
proportional to the pulse power and a second, integrated, signal 
proportional to the laser pulse energy. 

As shown in Figure 1, the power signal was used to trigger a 
Tektronix 7704 oscilloscope and was also displayed on that oscilloscope 
through the positive input in a differential amplifier. The energy 
signal was displayed on a Tektronix 564 b storage oscilloscope. 

The photodiode energy signal was calibrated by L. L. McKee [Ref. 38] 
using a Westinghouse RN-1 Laser Radiometer. The radiometer was used to 
provide an absolute measure of the incident pulse energy. 

The energy output of the laser was monitored on every shot and the 
power was displayed on the same oscilloscope used to take probe response 
data. For successive shots during the mapping procedure the energy was 
allowed to vary by + 10 %. Shots with energy variations exceeding this 
were excluded from the data base. When properly aligned the laser 
produced very consistent pulses. It should be noted that if shots were 



(5 Joules) power density of 2.8 X 10^ watts/cm 2 . 
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repeated at intervals of less than one minute, the laser energy tended 
to increase slightly from shot to shot. A one minute cooling period 
between shots resulted in excellent reproducibility. 

C. TEMPORAL REFERENCE FRAME 

The temporal relation between the diagnostic probe response and the 
leading edge of the laser pulse was determined by taking into account 
the various cable delays and laser beam optical path lengths. The 
geometry is depicted in Figure 1. A long (10 meter) delay line was 
used from the diagnostic probe to the oscilloscope in order to separate 
the laser pulse input and early diagnostic responses which would other- 
wise be superimposed by the differential amplifier. The computed delay 
was 33 nsec assuming the signal to travel at 1.98 X 10^0 cm/sec in the 
RG-17VU coaxial cable. 

To experimentally verify this delay, the 100 volt calibration signal 
from a Tektronix 555 Dual Beam Oscilloscope was fed into a Tektronix 
Tunnel Diode Pulser. The pulser output was fed into a "T." One side 
of the T went to the positive side of a Tektronix 7A-13 differential 
vertical plug-in on a Tektronix 7704 oscilloscope through the photo- 
diode power signal cable. The other side of the T went to the negative 
side of the same amplifier through the delay cable. A reproduction of 
the oscilloscopic trace obtained in this manner is given in Figure 2. 




37 nsec delay 
— 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 

Figure 2. Ten rreter cable time delay measurement. The horizontal scale 
is 20 nsec per division. The measured delay is 37 nsec. 
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When the measured 37 nsec delay is corrected for the difference in 
optical path length from the beam splitter to the target and the MgO 
diffuser block the delay is found to be 35 nsec. This experimentally 
measured delay agrees very well with that calculated earlier and is 
the delay used throughout this thesis in correlating diagnostic probe 
response with the laser pulse. Zero time is arbitrarily established 
as that time at which the leading edge of the laser pulse strikes the 
target. This is depicted in Figure 3- 



Figure 3- Definition of Zero Time. The horizontal scale is 20 nsec 



D. PLASM CHAMBER 

The target was located inside a chamber specifically designed by 
McKee [Ref. 38] to facilitate diagnosis of the laser plasma by probes or 
optical ireans. A horizontal cross section of this chamber is depicted 
in Figure 4. 




^ — Laser Pulse 
Defined Zero Time 



per division. 
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Lens 



/ 




Chamber Height 
is Six Inches 



Target Rotation Knob i» 



i 



Figure 4. Top view of vacuum chamber. Port #1 is the laser beam entry 
port, #2 is the reflected laser beam observation port, #3 is 
the transmitted laser beam observation port, ports numbered 
4, 5j 6, and 7 are optical/prcbe observation ports. There is 
also an optical/probe port in the top of the chamber. 
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The laser pulse enters the chamber after passing through a 
converging lens and strikes the target at an angle of 30 degrees with 
the target normal. This thirty-degree incident angle permits probing 
along the target normal which appears to be the axis of greatest 
symiBtry in relation to the expanding plasma. The converging lens was 
normally de focused to minimize cratering of the target. 

The target is a flat disk of one-half inch 6061 aluminum alloy 
plate. It is two inches in diameter. The target was rotated every 25 
to 35 shots to prevent severe cratering. 

The chamber was pumped to a vacuum of 2.5 X 10 3 torr (air) by an 
oil diffusion pump. This pressure was measured with an ionization 
gauge attached directly to the chamber. Unless otherwise stated, all 
measurements were made at this chamber pressure. 

^ Angular Reference Line 



Figure 5. Cylindrical polar coordinate system used in this thesis. 





Target Material 



Laser Focal 
Spot 
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The chamber coordinate system is polar with its center at the focal 
spot of the laser pulse. This is illustrated in Figure 5 - Throughout 
the text reference positions refer to cylindrical polar coordinates. 
Distances will be in centimeters unless otherwise specified. As an 
example, ( 0 . 5 , n /2 i 1 . 0 ) refers to the position r = 0.5 cm, 0 = n /2 and 
z = 1.0 cm. 



E. MAGNETIC PROBES 

Probe construction and calibration is covered by McLaughlin [Ref. 

393 * During this work a single magnetic probe consisting of five turns 

of Fbrmvar copper wire with a diameter of one millimeter was used. 

Using the procedures outlined by McLaughlin, an effective area of 
—6 2 

nA = 3*1 X 10 m was measured for the coil. The coil proved to have 
a linear frequency response up to 12 MHz. 

Since the probe was used only to determine temporal relations and 
to verify azimuthal symmetry of the self-induced magnetic fields, no 
integrator was used and the data was recorded simply as the voltage 
induced in the coil. 

F. FLOATING DOUBLE PROBES 

Appendix A contains details of the double probe construction and 
theory. The primary use of these probes was to determine the time at 
which disturbances arrived at the probe and to determine the relative 
density of the plasma at a given time in a given position. Probes were 
biased sufficiently negative to collect saturation ion currents. 

Probes used in this investigation have been numbered one through 
four. Probe number three has larger electrodes and was used 
extensively due to its higher signal to noise ratio. The characteristic 



19 



for the peak density of the main plasma pulse is reproduced as Figure 6. 
The characteristic was measured in a background of 2.5 X 10 J torr with 
a nominal laser energy of 5 +0.5 joules. The characteristic is 
symmetric and appears to partially saturate at + 10 volts. However, 
there is another, smaller, break in the curve at 20 to 25 volts. 
Therefore, the probe was biased at -20 volts to collect a saturated ion 
current. A large secondary emission coefficient of ( = 0.5 was used in 
obtaining a quantitative idea of plasma density to compensate for this 
large biasing voltage. 

By carefully isolating the entire system from ground, recorded 

probe response voltages at zero bias were reduced to less than 300 mV. 

For most probes the response was zero at zero bias indicating that the 

probe was in fact floating with the plasma potential. 

To obtain an estimate of the plasma density, equation (29) of 

Appendix A is solved for n^ and evaluated with Ze = 1.5j f = 0.5, 

R = 3 . 8 Q, A = 4.6 X 10“3 cm 2 and V given a variety of significant 

LP 

stream speeds: 

(1) n. (cm -3 ) = 1.59 X 10 20 V (volts) 

1 V (cm/sec ) 

J-UL 
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* 





Figure 6. Probe Nuinber 3 Characteristic at the Peak of the Main Plasma 
Pulse. 
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Table (1) gives versus probe response voltage for a number of 
significant plasma stream speeds. 

_3 

V LP (cm/sec) n_^ (cm )/Rrobe Response (volts) 

1.0 X 10 7 1.59 X 10 13 

5-0 X 10 7 3.18 X 10 12 

1.0 X 10 8 9-75 X 10 11 

Table (1) Plasma Density Versus Probe Response Voltage for Probe 
Number Three 

G. PRESSURE PROBES 

Pressure probes are useful in measuring the momentum per unit area 
transported by the plasma. In this investigation use was made of the 
piezoelectric properties of an X-cut quartz crystal. The momentum 
transfer caused by plasma impact with such a crystal creates a stress 
wave. Propagation of this stress wave through the crystal results in 
an induced voltage. Theoretical details describing this induced 
voltage are included in Appendix B. 

The pressure probe used in this investigation was constructed and 
provided by the Stanford Research Institute. 

H. DATA COLLECTION ERROR ESTIMATES 

1. Probe Positioning Errors 

The diagnostic probes were inserted into a brass collet system 
which was welded onto a brass base plate. The base plate could be 
moved about on an "0" ring seal located on the top of the chamber with- 
out breaking the vacuum seal. The base plate was scribed with 
reference lines. These lines allowed positioning of the plate over a 
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one millimeter grid attached to the top of the chamber. A new origin 
was determined at the start of each data session and marked on the grid. 
The probe was then moved to new positions using the grid as a reference. 
Hie position of the probe relative to the target was measured internally 
at the beginning of each data session and again at the end of the 
session to insure that the probe had not moved in its holder. While not 
very sophisticated, this method did allow positioning of the probe to 
within + 1/2 mm of the stated position in relation to the target . 

2. Absolute Plasma Density Errors 

The relative measure of plasma density from point to point is 
considered good due to the high reproducibility of the laser and plasma 
pulses. It is estimated that these relative values are accurate to with- 
in + 10 percent. 

Absolute plasma density calculations are subject to several 

errors . 

a. High Probe Bias 

A double probe bias higher than that at which initial 
saturation occurs may have generated errors as large as 100 percent in 
the recorded probe response due to secondary emissions and plasma 
perturbations. This error is consistent for all probe responses and 
can be ignored in the relative density plots. To help reduce the 
possibility of introducing a significant error in the absolute density 
determinations, a large secondary emission coefficient of (= 0.5 was 
used. 
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b. Deviations In the Laser Pulse Energy 

While the laser was very reliable and gave very reproducible 
pulses, the output energy was allowed to vary by as much as + 10 percent. 
The early signal probe response to the laser's energy content was found 
to vary linearly. The slope of the line describing this relation is 
such that the allowed deviation in laser energy could result in a 
deviation of + 7 percent in the amplitude of the early signal response. 

c. Error in Probe Tip Measurement 

The area of the probe tips is considered accurate to within 

ten percent. 

d. The plasma stream speeds are considered accurate to within 
ten percent. 

In view of these many deviations plus subtle factors [Fief. 48] 
in interpretation of double probe signals the absolute plasma densities 
are not considered accurate to better than a factor of two. 

3- Errors In Temporal Resolution 

The time resolution in these experiments is considered very good. 
At this background density, the plasma expansion speed is nearly 
constant within the allowed range of energies. The temporal resolution 
on individual photographs is enhanced by inclusion of the laser pulse 
which serves as a positive indication of zero tine. Temporal 
correlations and measurements are made with a horizontal scale of less 
than 40 nsec per cm. This allows interpretation of the results to with- 
in 4 nsec. For two double probes 1.75 cm apart, 4 nsec represents 46 

O 

percent of the transit time for the fast pulses travelling at 2 X 10° 
cm/sec and represents only 4 percent of the time involved in passage of 
the main plasma front at 2 X 10? cm/sec. When probe spacing and laser 



24 



energy errors are considered, the fast signal speeds should be good to 
within 47 percent and the main signal speeds to within six percent for 
single data points. In each case the speed has been averaged from more 
than 25 data points so that the respective limits are + 10 percent for 
the early signal speed and + 2 percent for the main plasma speed. 
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III. INDIVIDUAL EXPERIMENTS AND THEIR INTERPRETATION 



The thesis problem was attacked through a series of experiments 
designed to reveal the behavior and make-up of the early disturbance. 

A. DISTURBANCE SPEED MEASUREMENTS 

The first set of experiments was made in an effort to accurately 
determine the speed of component parts of the laser produced plasma. 
These speed measurements were accomplished in two different ways. 

1. Pulse Arrival Time at a Geographical Position 

In this experiment, a series of oscilloscopic photographs were 
taken at points along the z-axis. The series starts at 0.5 cm and 
consists of data points at one millimeter intervals to 3-0 cm, at two 
millimeter intervals from 3*0 to 4.0 centimeters and a final data point 
at 5 cm. This series was repeated three tines. 

a. The first series was taken with the oscilloscope set at 100 
nsec and 2 volts per cm. These settings are designed to enhance 
determination of the speed of the main plasma pulse. A typical probe 
response is traced as Figure 7 with the pertinent features labeled. 

Main plasma front passage times were measured for all data points. 

These tines are plotted against the probe position in Figure 8. The 
slope of the line connecting the data points represents the speed of the 
main plasma front at any point. As can be seen the velocity is nearly 
consistent at 1.1 X 10 7 cm/sec within the limits of this survey. 
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Figure 7- Double Probe Response at (0, 0, 1.0). The dashed line 

indicates the time at which the main plasma front passed 
the probe position. 




Figure 8. Time of Passage of the Main Plasma Front Versus Probe 
Position. 
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Assuming a constant speed and extrapolating backward Indicates that 
the main pulse left the target very early after laser pulse impact 
(approximately 10 nsec). 

b. The second series was taken with oscilloscope scales of 
20 nsec and 500 mV per cm. This scale enhances direct measurement of 
the early disturbances. A typical oscilloscopic trace is copied as 
Figure 9. 




Figure 9- Early Laser Plasma Produced Disturbances Recorded by a 
Double Probe at (0, 0, 1.25). 

The exact front of the early signal is often difficult to fix. 

This results in the large data scatter seen in early disturbance data. 

A plot of early disturbance front passage versus probe position is 
reproduced as Figure 10. The slope of this plot indicates a fairly 

g 

constant speed of 1.1 X 10 cm/sec. Also included in Figure 10 is a 

similar plot for the second early disturbance. The slope of this line 

7 

indicates a front speed of 5*9 X 10 cm/sec. 

It should be noted that the speed of these early disturbances is 
reasonably constant even to the point where the signals are lost in 
system noise. 

By extrapolating backward assuming a constant speed we find that 
the various disturbances left the target at the following times; main 
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plasma at t - 0 to 10 nsec, first early disturbance at about 3 nsec and 
the second early disturbance at t = 12 nsec. 




Figure 10. Plot of the Time of Passage of the Early Disturbance Versus 
Probe Position 

2. Two Double Probes 

To obtain more accurate speed information, a two probe holder 
was constructed to allow simultaneous sampling of the same plasma by 
two tandem probes. The holder is illustrated in Figure 11. 

The probes used were very fine to minimize perturbation of the 
streaming plasma. The probes were placed in the holder 1.75 cm apart. 
They were positioned in the chamber so that one probe was one millimeter 
to the right of the target normal and the other was one millimeter to 
the left of the target normal. Each probe was biased at -20 volts and 

_5 

the chamber pressure was reduced to 1.8 X 10 torn. A series of 32 
pictures were taken on a Tektronix 555 Eual Beam Oscilloscope. The 




Insertion Collet System Movable Brass Plate 





Figure 11. Device for the Insertion of Two Double Probes 

scope was set at 500 mV and 40 nsec per cm. A typical oscilloscopic 
trace of the probe response is depicted in Figure 12. At ten of the 
data points, the probe cables were reversed. In each instance, the 
same separation between arrival times was noted. The reversed cable 
trace is also included in Figure 12. 
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probe l 

probe 2 



Figure 12. Simultaneously recorded tracings from two double probes 
sampling the same plasma, (a) The first probe is at 
(0.1, tt / 2 , 1.1) and the second probe is at (0.1, 3^/2, 

2.85). (b) Taken at the same data points but with the 

cables reversed at the oscilloscope. 



The average speed of the early signal as computed from this data 

g 

base is 1.0 X 10 cm/sec. The average speed of the main plasma signal 

7 

is found to be 1.1 X 10 cm/sec. These results agree very well with 
those found by the first method. A summary of these speeds is given 
in Table 2. 



Plasma Feature 


Speed 


First Early Disturbance 


1.1 X 10 8 cm/sec 


Second Early Disturbance 


5.9 X l 0 7 cm/sec 


Main Plasma Front 


1.1 X l 0 7 cm/sec 



Table 2. Summary of Significant Speeds 
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B. EFFECT OF A TRANSVERSE ELECTRIC FIELD 



A variable transverse electric field of up to 13,000 volts/cm was 
applied to a collimated plasna to determine whether the field would 
deflect the plasma beam. This will reveal information on the charge 
to mass ratio of the disturbance if it is particulate in nature. 

Figure 13 depicts the device constructed to provide a collimating 
effect on the laser produced plasma. The devise was designed with 
several purposes in mind: 

1. To collimate the plasma. 

2. To help shield the electrodes from the plasma and thereby 
reduce the possibility of discharge across the conductive plasma. 

3. To reduce the total number of plasma particles subject to the 
field in order to reduce the effects of shielding the plasma. 




Figure 13. Collimating Device (a) General View, (b) Top View, and 
(c) Front View 
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Figure 14 depicts the probe response at (0, 0, 5*8) with and without 
the plexiglass collimating device in place. As can be seen, the probe 



response has been altered; however the early disturbance and main plasma 
pulses remain distinguishable features of the response and therefore the 



pinhole collimating device is usable. 




100 nsec/soomV 



Figure 14. Effect on the plasma produced by the presence of the 

collimating device, (a) Double probe response without the 
device, (b) Double probe response with the device in place. 
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Figure 15. Electrodes and Their Mounting 
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Figure 15 depicts the electrodes used in this experiment and the 
mounting equipment. All surfaces of the entire structure are heavily 
coated with General Electric "RTV" to eliminate corona formation. 

Assuming that the electrodes form a pair of infinite parallel 
plates (ignore fringing), the stored charge on these electrodes is found 
to be: 

(2) q = EA Where (A) is the plate area 

(3) V = Ed Where (d) is the electrode spacing 

(4) q = VA/d 

2 

With V/d = 13,000 volts/cm and (A) = one cm the stored charge is 
1.5 X 10 ^ coulombs or 7-2 X 10^ electrons. 

Plasma densities greater than this will be able to completely shield 
themselves. Typical electron densities in the main plasma pulse are 

i/j _3 

approximately 10 cm while probe response to the early disturbance 

indicates a density of approximately 2.4 X 10 12 ion charges/cm^. 

Assuming a nearly neutral plasma in each case it is obvious that both 

pulses are capable of completely shielding themselves. The collimating 

device has pinholes 3-2 millimeters in diameter. This reduces the early 

11 -3 

pulse charge density to 2.5 X 10 cm . The early plasma is still 
capable of shielding itself. However, in doing so, nearly all of the 
particles would be lost from the beam resulting in a significantly 
reduced signal. This argument assumes that the plasma particles have 
time to respond to the applied field. 

Because of the shielding problem a large capacitor was used in the 
electrode charging circuit. This circuit is shown in Figure 1 6. 
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Figure 16. Electrode Charging Circuit 



The RC time for this circuit is approximately 2 nsec between the 
storage capacitor and the electrodes. The capacitor stores 4.5 X 10“^ 
coulombs of charge at 18,000 volts. Therefore the electrodes will 
receive a significant amount of additional charge from the power supply 
during the 27 nsec in which the early disturbance is traversing the field. 
The charge available in the storage capacitor should be sufficient to 
overcome any shielding in the early signal. 

Deflection of singly ionized aluminum atoms assuming no shielding 
of the field by the early disturbance will be: 

(5) D = at 2 /2 Where (D) is the lateral displacement 

(6) a = F/m = eE/m = eV/md 

If V/d = 13,000 volts/cm then the acceleration is 4.6 X lO 1 ^ cm/sec 2 

17 2 

for aluminum ions and 8.5 X 10 cm/sec for electrons. 

The early disturbance travels with a speed of approximately 
1.1 X 10^ cm/sec. Therefore, it will be within the 3 cm field for 
t = 27 nsec. Substituting these values into equation (5) we find that 
the ions incur a transverse displacement of 0.9 mm during transit 
while the electron lateral displacement is 170 cm (ignoring space charge 
effects) . 

These values indicate that if the early disturbance was cortposed 
entirely of electrons, it should be completely deflected and no signal 
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would be received. However, if the early disturbance Included ions, 
then these ions would not achieve sufficient deflection to miss the 
exit pinhole and any signal attenuation would be very small. 

An unforeseen problem arose in this experiment either because of 
the large distance at which the probe samples the plasma or because of 
the pinhole collimator. The plasma sampled by the probe was not 
reproducible. Major features were always identifiable but their 
amplitudes differed by as much as a factor of three from shot to shot. 
This unresolved problem has greatly complicated the analysis by causing 
a large deviation in the data. 

Results of this experiment are reproduced in Figure 17 . The straight 
line is a least squares fit to the data. Its slope is 0.19 mV/kV 
indicating that there is little dependence of the early disturbance on 
the applied transverse field. These results suggest that the early 
signal is not composed entirely of electrons but that there are heavier 
particles present. 
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Figure 17 . Early Disturbance Peak Amplitude Versus Transverse Electric 
Field Strength. 
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C. EFFECT OF A TRANSVERSE MAGNETIC FIELD 



Because of the complications inherent in the production and inter- 
pretation of data taken with a large transverse electric field, a 
transverse magnetic field was applied to the streaming disturbances. 
The experimental arrangement is depicted in Figure 18. 




Figure 18 . Experimental setup for the application of a transverse 
magnetic field. 



The magnetic field was measured with a gaussmeter. The field in 
the center of the magnetic pole area is 950 gauss. Near the surface of 
the poles the field strength increases to 1325 gauss. 

Great care was taken to insure that the incoming laser pulse did 

n 

not strike the target and that the distance between the magnetic poles 
was sufficient to allow plasma expansion without perturbation. In 
order to determine if the magnetic field was inducing currents in the 
probe a number of pictures were taken with the target removed allowing 
the laser pulse to transit the chamber and exit through port number 
three. Figure 19 shows the oscilloscope response in this situation. 
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Figure 19. Double probe response in the presence of a 950 gauss 
magnetic field. The target has been removed and the 
laser pulse transits the chamber and exits through a 
glass cover over port number three. 

These tracings indicate that no probe response is generated by either 
the magnetic field or laser pulse interaction with components of the 
apparatus other than the target. 

Gyro frequencies for electrons and ions are given by equation (7): 

. (7) oj = eB/mc (radians/sec) 
c 

The field in this case is 950 gauss so that the ion gyrofrequency 

C in 

is 3*35 X 10^ rad/sec and the electron gyrofrequency is I .67 X 10 
rad/sec . 

The Larmor radii are given by equation (8) 

( 8 ) = v ± /„ c 

From this expression we find that for electrons the Larmor radius is 
9.75 X 10“3 cm and for aluminum ions it is 48*J cm. From this it is 
obvious that if the early disturbance is composed essentially of 
electrons, the signal should be greatly diminished at distances greater 
than a few gyro radii or 0.5 cm. However, if the early disturbance is 
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essentially a neutral plasma then there should be little change in the 
amplitude. The ions with their much larger gyro radius will hardly be 
deflected by the field. As soon as sufficient space charge separation 
is created, the ions will drag the electrons out of the magnetic field. 

Other researches [Refs. 1, 6, 25, 28 and 33] have reported that 
plasma fronts tend to be rich in electrons while the main body is 
ionically rich. In this case, the front will be changed from electron 
rich to ion rich by the magnetic field. If we consider the plasma on 
a microscopic, interparticle, basis and set the Lorentz force equal to 
zero we can obtain a measure of the space charge separation necessary 
to pull the electrons out of the magnetic field 



(9) 


V^p X "b = E 


= VppB since ^p_L "B 


(10) 


Vpp B = ke/r 2 


which on solving for r gives : 


(ID 


1/2 

r = (ke/VB) 


where k = 1 



If we substitute V TD = 1.1 X 10 7 cm/sec and B = 950 gauss into equation 

i-LL 

(11) we find that the separation is 0.1 mm. 

Figure 20 duplicates photographs of the oscilloscopic record of the 
double probe response with and without the transverse magnetic field. 
There is no discernible change in the early signal when the magnetic 



field is introduced. 



(b) without magnetic field 

< a > with magnetic field 




Figure 20. Double probe response at (0, 0, 1.25) to the streaming 
plasma with (a) and without (b) a transverse 950 gauss 
magnetic field. 

It is interesting to observe that the main plasma is delayed by 
about 10 nsec and the early plasma by about one or two nsec. For the 
main plasma traveling at 10^ cm/sec this delay represents a change in 
plasma front position of one millimeter. This compares with the 
prediction obtained from equation (11) if the electron rich area 
originally extends 0.9 mm in front of the neutral boundary of the plasma. 

It might be argued that since the direction of this transverse 
magnetic field is horizontal there is a drift velocity due to 
the gravitational force. Depending on the direction of the fields 
this drift could be either back into the target or away from the target 
toward the probe. The magnitude of this drift is given by: 
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U2) V ra 



me g X B = meg/qB 
2 



qB 



Fbr the values under consideration in this experiment the magnitude 
of VpQ is 5.9 X 10-8 cm / sec which is insignificant. 

Since the magnetic field has little effect on the early disturbance 
this experiment strongly suggests that the early disturbance is not 
purely electronic in composition. If it were, then the magnetic field 
would deflect the electrons into a gyro radius with diameter much 
smaller than the distances at which the probe measurements were made. 
This would result in the loss of the early disturbance and the probe's 
response would drop significantly. Therefore, if the early disturbance 
is particulate in nature, it must contain heavier ions. 



D. SEPARATION OF THE PHOTOELECTRIC RESPONSE 

During the density contour mapping a disturbance of small half width 
(4 nsec) was observed to occur at approximately t = 16 nsec, i.e. near 
the peak of the laser pulse. The time of occurrence was constant for 
all points on the mapping grid indicating a very high propagation speed. 
Several authors [Refs. 16 and 33] have reported strong X-ray and UV 
radiation emanating from the target area of a laser produced plasma. 
Ready [Ref. 40] has ascribed the early disturbance noted on his 
Langmuir probes to reverse photoelectric effect. His measurements were 
made at 30 cm from the target in near vacuum. As will be noted later, 
the early plasma decays rapidly between 0.75 and 2.0 cm. Therefore at 
30 cm it is most probable that all he saw in terms of early probe 
response was the reverse photoelectric response. At distances less 
than one centimeter from the target the photoelectric response is masked 



42 



by the much larger response due to the early plasma pulse. Figure 21 
Illustrates this masking and shows the appearance of the photoelectric 
signal at 1.7 and 3.0 cm. 




Figure 21. Double probe response to a streaming laser produced plasma 
at (a) (0, 0, 1.3), (b) (0, 0, 1.7) and (c) at (0, 0, 3-0). 
Asterisks mark the suspected photoelectric response 

In order to verify that this signal was photoelectric in origin, a 
fused quartz (SiO^) plate was obtained. The plate is 2.5*1 cm square and 
and has a thickness of 0.02 cm. The transmissivity of this plate is 
greater than ten percent for wavelengths longer than 0.12 microns. 

Figure 22 depicts the transmissivity as a function of wavelength for a 
one-half inch plate of the same material. The work function for tungsten 
is approximately 4.55 eV. Therefore, the maximum wavelength for single 
photon photoelectric emission from the probe tips is 0.267 microns. This 
wavelength is marked on Figure 22. The hashed area denotes the trans- 
mitted wavelengths capable of producing the observed photoionization. 



IOOt 




Figure 22. Transmissivity of a 1/2 inch fused quartz plate. The 
hashed area represents those wavelengths that are 
transmitted and are capable of ionizing tungsten.' 



In this experiment, the quartz plate was placed between the probe 
and the target. Probe response with and without the quartz was recorded 
and a comparison made. These oscilloscopic traces are reproduced in 
Figure 23. 




Figure 23- Oscilloscopic record of a double probe response at (0, 0, 
1.75) with (a) the laser beam passing through the chamber 
but no target in place, (b) without the quartz plate 
screen and (c) with the quartz plate screen in place 
between the target and the probe. An asterisk marks the 
verified photoelectric response 
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This experiment was repeated ten times at different locations in 
the chamber with identical results. As seal, the quartz plate stops 
all physical particles but passes UV radiation which produces the 
photoelectric response at the probe. 

The height of the response is approximately 300 mV. Assuming that 
ten percent of the electrons ejected from the negative electrode are 
collected at the positive electrode this corresponds to a current flow 
of 0.079 amps for 4 nsec or a total of 1.97 X 10^ electrons. Assuming 
an average photon energy of six eV and also assuming that ten percent 
of the incident photons result in the ejection of an electron from the 

negative electrode we have an energy deposition on the probe tip of 

—8 

1.9 X 10 joules. The cross sectional area of this tip is 
_o 2 

4.6 X 10 D cm . Assuming that the electromagnetic radiation is trans- 
mitted equally in all directions through a hemisphere whose center is 
the focal spot on the plane target, this corresponds to a total energy 
of 7.9 X lO - '* joules transmitted at wavelengths between 0.12 and 0.27 
microns. This represents about 1.3 X 10 of the incident energy in the 
laser pulse. 

Equation (25) predicts that for a static plasma the ion density 
is [Ref. 45]: 

(25) n. = 1.67 I + (m 1 /kT e ) 1/2 /eA p (em“ 3 ) 

If the 300 mV probe response is due to photo ionization of the background 

this would correspond to a background gas ion density of 

1.79 X lO 1 ^ cm _ 3 assuming kT e = 3 eV, rr^ = 28.96 AMU and that the total 

-2 2 

surface area of the probe electrode is 1.4 X 10 cm . Since the 
background gas density is only 8.85 X 10^ cm ^ this response must be 
due to the reverse photoelectric effect described by Ready . 
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E. DENSITY MAPPING OP THE EARLY PLASMA 



The expanding plasma density has been mapped using a double probe. 
The mapping grid is shown in Figure 24. 




Z ► 



Figure 24. Density Mapping Grid 

The general method used was to acquire an oscilloscopic record of 
the probe response at each grid point. A zero time was marked on each 
photograph and corrected for the 35 nsec cable delay . Probe response 
voltages were then read from the photograph corresponding to the 
desired time intervals for the contour maps. These voltages are readable 
to within + 0.2 volts. The voltage from each grid point corresponding 
to a specific time is then plotted on the grid depicted in Figure 24. 
Contours were drawn by hand. All of the data and contour mappings were 
analyzed manually. However, because these are only relative density 
contour maps and because of the many data points involved in each map, 
any errors introduced should be minimal. 
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The area excluded from each map is the corridor through which the 
laser pulse travels in reaching the target. Since plasma density is 
proportional to probe response voltage, these contour maps give 
relative density contours. 

For the early plasma pulse, the probe response was mapped at 10 nsec, 
20 nsec, 30 nsec, 40 nsec and then at 20 nsec intervals to 160 nsec. 

This allows a detailed analysis of the expanding plasma pulse. These 
contour maps are included as Figures 26 to 35. 

The maps indicate that the early disturbance expands rapidly 

Q 

(10 cm/sec) along a line inclined at approximately 30 degrees with the 
target normal. This line also describes the path followed by the 
reflected laser pulse. 

It should be noted that typical densities measured throughout these 

12 -2 

maps are on the order of 10 cm . This is nine orders of magnitude 
below the critical density so that the plasma being examined here is 
essentially transparent to the incoming laser radiation. 

The focal spot size of this laser pulse is 1.5 mm in diameter. It 
should be noted that the early disturbance expands very rapidly to a 
full radius of approximately 1.0 cm by z = 0.2 cm. At z = 0.2 cm the 
plasma radial expansion stops and all further expansion occurs along 
the axis of symmetry. This behavior is strongly suggestive of a force ^ 
acting on the plasma to pinch the radial development. Since no 
external fields other than the small probe potential were applied, it 
is evident that this force must be produced internally within the plasma. 
The self-induced magnetic fields mapped by McKee and Bird [Refs. 8 and 38] 
would provide a pinching effect on the plasma. The J X B force was 
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neglected by McKee as too small to be significant. However, further 
investigation of this force seems warranted and is included in Section 

IV. B. 

It should also be noted that the plasma density gradients are 
higher in the radial direction than in the direction of pulse 
propagation. This is also strongly suggestive of a radial confining 
force. These gradients are plotted in Figure 25. 

Maximum densities noted in the early pulse are approximately 
4 X 10 cm . This corresponds to the 2.5 volt contour in each of 
Figures 26 to 29. In Figure 30 the early plasma pulse has separated 
from the main plasma pulse. Subsequent figures contain only main plasma 
contours. The maximum early plasma densities are greater than the 
background gas density by a factor of five. 

It should be noted that as early as t = 10 nsec there is evidence 
of plasma at distances as great as 0.8 cm. This plasma may partially 
be the result of intense background gas photoionization. However, the 
magnitude of the response indicates ion densities three orders of 
magnitude higher than the background gas density (assuming a 
macroscopically quiescent plasma) . 

These topographical density maps also show that after flooding into 
the chamber, this early disturbance decays rapidly so that at 100 nsec 
there is little remaining evidence of the early plasma pulse. Another 
interesting feature is that the decay appears to proceed more slowly 
along the target normal than along the axis of propagation of the early 
disturbance (see Figures 30 and 31)* 
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Early plasma pulse density gradients in (a) the direction 
of propagation and (b) in the radial direction 
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Figure 26. Laser-produced plasma relative density contours at 
t = 10 nsec. 
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Figure 27 • Laser-prcduced plasma relative density contours at 
t = 20 nsec. 
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Figure 28. Laser-produced plasma relative density contours at 
t = 30 nsec. 
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Figure 29. Laser-produced plasma relative density contours at 
t = 40 nsec. 
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Figure 30. Laser-produced plasma relative density contours at 
t = 60 nsec. 
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Figure 31* laser-produced plasma relative density contours at 
t = 80 nsec. 
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Figure 32. laser-produced plasma relative density contours at 
t = 100 nsec. 
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Figure 33* Laser-produced, plasma relative density contours at 
t - 120 nsec. 
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Figure 3^. Laser-produced plasma relative density contours at 
t = 140 nsec. 
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Figure 35. Laser-produced plasma relative density contours at 
t = 160 nsec. 
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F. DENSITY MAPPING OF THE MAIN PLASMA 

As an extension of section IV. E. , relative density contour maps 
were completed at 100 nsec intervals from 100 to 800 nsec. The sane 
mapping grid as shown in Figure 24 was used to produce these maps. 
However, the oscilloscope scales were increased to 100 nsec and 2 volts 
per division. The results are illustrated in Figures 37 to 44. 

The first, very obvious feature of the main plasma is noted at 60 
nsec in Figure 36. The main plasma body with densities corresponding 
to probe response voltages greater than 14 volts has been shaded. This 
map indicates that at early times the main plasma does not expand 
symmetrically about the target normal. The radial density is strongly 
skewed to the left between z = 0.2 and z = 0.8 cm. The high density 
area appears to be centered at (0.6, tt/ 2, 0.5). Since these contour 
maps were made for only a single plane, no determination of the three 
dimensional distribution is possible. Examination of Figures 30 to 34 
in section IV. E. indicates that this asymmetry decays rapidly so that 
by 120 nsec the main plasma has become essentially symmetric about the 
target normal. To further investigate and hopefully verify the late 
time symmetry of the main plasma with respect to the target normal, a 
special target of 1/8 inch thick 6061 aluminum alloy was constructed 
and mounted onto the face of the half inch target normally used. A 
28 cm converging lens was used to focus the laser pulse to the smallest 
obtainable focal spot size at the target face. This was accomplished 
with an autocollimator. 

Two double probes (numbers one and two) were positioned inside the 
chamber as depicted in Figure 45. 
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Figure 36. Early time asymmetry of the main plasma with respect to 
the target normal 
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Figure 37* Laser-produced, plasma relative density contours at 
t = 100 nsec. 
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Laser-produced plasma relative density contours at 
t = 200 nsec. 
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Figure 39* 



Laser-produced plasma relative density contours at 
t = 300 nsec. 
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Figure 40. Laser-produced plasma relative density contours at 
t = 400 nsec. 
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t = 500 nsec. 
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Figure 42. Laser-produced plasma relative density contours at 
t = 600 nsec. 
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Figure 4 3- Laser-produced plasma relative density contours at 
t = 700 nsec. 
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Figure 45. Experimental arrangement for verification of late time 
plasma expansion symmetry with respect to the target 
normal. Probe one is at (0.2, n /2, 2.5) and probe two 
at (2.2, 3n/2, 2.5). 



A total of 50 shots were made without disturbing the target. 

Pictures of the double probe responses were taken every five shots as 
displayed on a Tektronix 555 dual beam oscilloscope. The laser had 
burned a visible hole into the target after the fourth shot. 
Reproductions of the oscilloscopic traces at ten shot intervals are 
included as Figure 46. 

Some increase in the signal at probe number two was expected due 
to a shotgun effect as the plasma expands from the base of the deepening 
hole. However, as the traces in Figure 46 indicate, there is no 
appreciable change in the plasma density sampled by either probe during 
the entire sequence except that the early plasma pulse is somewhat 
attenuated after the fifth shot. This behavior strongly suggests that 
the expansion dynamics are independent of the plasma behavior prior to 
its rapid expansion into the background. 
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Shot 1 



Shot 30 




Shot 10 



Shot 20 



All Scales 100 




Shot 40 



Ia) Shot 50 



Figure 46. Series of oscilloscopic traces taken on a dual beam 

oscilloscope during the burning of a hole 1.029 cm deep 
in an aluminum target . (a) probe two and (b) probe one 

t 

Following this experiment, a detailed examination of the crater was 
undertaken with an electron microscope. Figure 47 is an electron 
micrograph of the entry point of the laser pulse, Figure 48 depicts 
the exit hole from the 1/8 inch plate. Figure 49 shows a cross section 
of the hole in the 1/8 inch plate. Figure 50 is a cross section of the 
crater burned into the 1/2 inch back plate. Figure 51 shows a cross 
section of the base of the hole excavated into the back plate and 
Figure 52 depicts a typical crater wall. 

Careful measurement of the volume of the cavity reveals that 
approximately 0.014 cm^ of aluminum were excavated. This represents a 
total of 8.44 X 10^° aluminum atoms or approximately 1.7 X lO 1 ^ atoms 
per laser pulse. This value is two orders of magnitude larger than 
the 10 1 ? atoms reported by Honig [Ref. 23]. However, examination of 
Figure 47 indicates evidence of splashes of molten material radiating 
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Figure 47. Entry Hole in the 1/8 inch Aluminum Target- 
Magnification is 18X. 




Figure 48. Exit hole as the laser pulse leaves the 1/8 inch tar t 
plate. Magnification is 40X. 
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Figure 49 . 



Cross section of the hole burned through the 1/8 inch 
aluminum target pla'e. Magnification is 23X. 



Figure 50. 



Cross section of the crater burned into the 1/2 inch 
aluminum back plate. Magnification is 17X. 
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Figure 52. Typical section of the crater wall in the 1/2 inch aluminum 
back plate. Magnification is 540X. 



Figure 51. Base of the crater burned into the 1/2 inch aluminum tack 
plate. Magnification is 480X. 
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ft?oni the entry point of the laser pulse. Examination of the probes 
after the experiment reveals that each is heavily coated with a layer 
of shiny aluminum. These splashes and the aluminum coating have not 
been observed in previous experiments using the laser in a defocused 
mode. This evidence indicates that a great deal more material was 
ejected during the hole burning process than is ejected when the laser 
pulse is defocused to produce only a shallow crater. Rather than being 
ionized, this additional material appears to be essentially neutral and 
may even be in a partially condensed state. This statement is 
supported by the aluminum coating on probe number two which did not 
indicate the presence of a significant number of ions at any time during 
the experiment. 

An obvious cause for the increase in ejected material is that the 
hole tends to contain all of the incident laser energy whereas in the 
shallow crater produced by the defocused laser, a significant amount 
of the incident energy may be reflected from the target and expanding 
plasma. Examination of Figure 51 reveals evidence of dynamic shock 
processes in the base of the crater. Large pieces of aluminum are 
apparently being ejected by shock waves. The ragged appearance of the 
entire surface of the crater suggests that dynamic removal of material 
may be more important than thermal processes. Inhomogeneities in the 
aluminum alloy structure may strongly contribute to the appearance. 

This subject will not receive further consideration in this 
investigation. 

Main plasma density gradients in the radial and axial direction are 
plotted in Figure 53 for times of 100 nsec, 400 nsec and 800 nsec. The 
radial density gradient is plotted at z = 1.0 and z = 2.0 cm. The axial 
density gradients are along the target normal. 
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These main plasma density contours lend further evidence for radial 
confinement of the plasma. This evidence is seen in the cylindrical shape 
of the expanding plasma out to 3*0 cm (see Figure 40) and by the 
initially steep and slowly diffusing radial density gradients seen in 
Figure 53- 

The high density gradients for small z values at long time intervals 
seen in Figures 43 and 44 must be the result of rapid plasma cooling in 
contact with the target face since the laser pulse shuts off at 
t = 30 nsec. 

To obtain a more detailed picture of the plasma density gradients 
close to the target under slightly different conditions, data was collected 
on a one millimeter grid from z = 0.2 cm to z - 1.2 cm and radially to 
1.5 cm. Contour maps were constructed at 100 nsec intervals from 200 nsec 
to 800 nsec. These maps are included as Figures 55 and 56. The data were 
obtained with probe number three biased at -15 volts and with the laser 
energy reduced to 4.0 +0.4 J. 

The most striking feature of the plasma is that it retains nearly 
the same density gradients and radial dimension throughout the tine scale 
of this survey. This is most dramatically illustrated in Figure 54 where 
the one volt contour for 200 nsec is superimposed on the one volt contour 
for 800 nsec. As can be seen in Figure 54, the radial extent has 
increased from 0.7 cm to only 1.0 cm during this time interval. This 
provides further evidence for the presence of a radially confining force. 
Also note that between 0.5 cm and 1.2 cm the density gradients are almost 
entirely in the negative radial direction so that any temperature gradient 
in the negative axial direction would provide the necessary source term 
for the reversed self-induced magnetic field reported by Bird. 
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Radial Gradients 




Figure 53. Main plasma density gradients at (a) 100 nsec, (b) 400 nsec 
and (c) 800 nsec. 
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Figure 5^- Density contours for (a) 200 nsec and (b) for 800 nsec. 



78 





1 

Axial Distance (cm) ► 



Figure 55. Laser produced plasma density contours at (a) 200 nsec, (b) 
300 nsec, (c) *100 nsec, (d) 500 nsec and (e) 600 nsec. 
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Figure 56. Laser produced plasma density contours at (a) 700 nsec 
and (b) 800 nsec. 
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G. DEPENDENCE OF THE EARLY PLASMA ON LASER PULSE ENERGY 

The maximum density in the early signal has been found to vary 
linearly with the laser pulse energy. A series of oscilloscopic traces 
of the double probe response at (0, 0, 1.9) were recorded with the 
laser pulse energy varying from 2.9 joules to 9.2 joules. The peak 
early probe response versus laser energy is plotted in Figure 57* 




Figure 57. Peak Double Probe Response to the Early Plasma Versus Laser 
Pulse Energy 

It appears that the early response has a linear dependence on the 
laser energy. For this particular probe (probe number three) and 
within the energy range plotted, the dependence can be expressed in the 
relation : 

(13) V pp = (l.lE^ - 0.375) volts 
Where V pp is the peak probe response to the early plasma pulse and E^ 
is the laser pulse energy in joules. It should again be noted that these 
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data and equation (13) are valid only for this experiment. The 
amplitude of the early pulse is highly dependent on background gas 
pressure as shown by Arifov [Ref. 5]. These data were collected in 
near vacuum at 10 J torr. This pressure is at the threshold of the 
pressure dependent region described by Arifov. There are peak ion 
densities of 2.5 X 10^ cm in the early plasma pulse at this pressure. 
As the early plasma pulse propagates, it ionizes the background. 

However, if the background density is less than the ejected plasma 
density then the additional ionization will be comparatively small. 

The background density is equal to the early plasma density at a 

_5 

pressure of 7 X 10 torr. As the background pressure increases above 
this value, ionization by the streaming plasma increases the total ion 
density detected by the probe and the pressure dependent region is 

_5 

reached. The significant boundary at 7 X 10 torr is consistent with 
Arifov 's and Isenor's work [Refs. 5 and 25]. 

H. TIMING SEQUENCE FOR DISTURBANCES LEAVING THE TARGET 

An attempt has been made to determine the times at which the various 
disturbances left the target. This was accomplished by assuming that 
the speed of the disturbance is constant between the target and the data 
point. By noting the time at which the disturbance reaches a specific 
data point and subtracting the transit time using the measured speed of 
the disturbance an estimate of the time at which the disturbance left 
the target is obtained. This method appears to give reasonable results 
for the early disturbance and the photoelectric response due to the 
constancy of their propagation speeds. Repeated measurements with 
available data were averaged. The average computed departure time is 
t = 16 + 2 nsec for the photoelectric disturbance and t = 18 +6 nsec 
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for the early plasma disturbance. The average of the times at which 
the main plasma left the target is t = 2 + 15 nsec assuming a constant 
speed of 1.1 X 10^ cm/sec. This last value does not seem consistent 
with any dynamic model for the plasma generation. It is assumed that 
the error is introduced because the main plasma speeds are much higher 
close to the target. It is not possible to measure the position of 
the main plasma front at distances less than 0.6 cm because the early 
disturbance is still separating from the main plasma body. This is 
illustrated in Figure 58. This problem plus the higher oscilloscope 
time scale (100 nsec per division) used to measure the main plasma 
account for the high standard deviation in this computation. 




Figure 58. Double probe response at (0, 0, 0.6) showing incomplete 
separation of the early disturbance from the main plasma 
body. 

I. CORRELATION OF THE QUARTZ PRESSURE PROBE AND DOUBLE PROBE RESPONSES 
In this experiment a double probe and a quartz pressure probe were 
placed side by side at an axial distance of 1.4 cm from the target. 
Probe responses were recorded from each diagnostic probe on successive 
shots using the same delay cable. 
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In this investigation, the momentum transport of the plasma provides 
the strain on the quartz pressure probe. Assuming that plasma particle 
collisions with the pressure probe are elastic and that the mass of the 
plasma is much less (6 orders of magnitude at least) than that of the 
probe and its holder, the force density producing the strain in the 
gauge is: 

(13) F - dP - + Ve n e + VnV 

Since the electron mass is negligible compared with the atomic mass we 
have m^ = ir^ » ' m g and equation (13) can be approximated by: 

UH) P - d_(”V 1 V n i + "n) > 

It should be noted that this equation assumes that the plasma flow 
velocity, V^p, is much larger than the plasma ion thermal speed. This 
assumption is justified in Appendix A. 

Applying the assumption of elastic collisions and doing a simple 
dimensional anaylsis on equation (14) indicates that: 

(15) F = A P = 2mV Lp (n i + n^/T^ 

At 

Where T^p is the time for the plasma to travel one cm. Since 
Tpp = (Vpp) -1 equation (15) can be rewritten as: 

? 3 

(16) F = 2m(Vjp) (nj_ + r^) newtons/cm 

This analysis describes the plasma as a square pulse which it is 
not. However, for the purposes of this discussion the approximation 
of a square plasma pulse will do. For the main plasma pulse with 
m = 27 AMU, = 1.1 X 10 7 cm/sec and (n ± + r^) = 2 X 10 15 cm -3 , the 
strain producing force density is 0.08 newtons/cm 3 . 
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A similar analysis of the early plasma pulse with V TT3 = 1.1 X 10® 

Lll 

cm/sec and (n^ + r^) = 2.5 X 10 -*-^ cm“3 yields a strain producing force 
density of 0.01 newtons/cm 3 . 

The ratio of these two stress producing force densities is 0.125. 
This is also the ratio of voltage outputs from the crystal in the 
presence of each pulse since the probe response is proportional to the 
stress. 

These impulsive forces produce stress waves in the aluminum gauge. 
These waves propagate with the speed of sound. For aluminum, 

= 6,420 m/sec [Ref. 4]. Therefore it takes 390 nsec for the stress 
wave to propagate through the 0.25 cm aluminum cover disc and epoxy to 
the quartz gauge. 

Figure 59 depicts the signals produced from the same position by 
(a) the double probe and (b) the pressure probe. As can be seen, the 
main plasma front corresponds very well with the large, late response 
from the pressure probe. However, the early response from the pressure 
probe occurs approximately 200 nsec before a stress wave generated by 
the early plasma would arrive at the crystal. No conclusive explanation 
is offered for this temporal discrepancy in the arrival of the early 
pressure probe response. Perhaps the speed and sharpness of the early 
plasma cause the wave to travel through the aluminum at a speed in 
excess of the speed of sound. The electrodes are well shielded with 
epoxy and the outer case is isolated from the circuit so that electro- 
static pickup is unlikely. However, the early signal does correlate 
with the time an expected electrostatic pickup would be measured. It 
is also possible that the obvious early signal is not that associated 
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with the early momentum transfer. The expected amplitude of the early 
signal is approximately 4 mV. A signal of this magnitude could easily 
be lost in the noise. 

The fact that there are two distinct disturbances registered by the 
pressure probe indicates the presence of two momentum transporting pulses. 
However, the poor temporal correlation of the early signal makes the 
results of this experiment somewhat inconclusive. 

It should be noted that researchers at the Stanford Research 
Institute have fired high power, 30 nsec, laser pulses at the gauge's 
aluminum target face. The quartz gauge then measures the momentum 
transported to the target by the fast blow off of laser produced plasma. 
These investigations have revealed two momentum impulses delivered to 
the target as depicted in Figure 60. 




100 nsec/500 mV per div. 



pressure probe response 
target removed 



100 nsec/20 mV per div. 



time main plasma impulse expected 
time early impulse expected 

early pressure probe response 



Figure 59. Correlation of signals received by (a) the double probe 

and (b) the pressure probe at (0, 0, 1.4) in the presence 
of a laser plasma. 



86 



4 




fso laser impacts Al Cover 



Figure 60. Stress time history In a laser Irradiated sample in 
vacuum. r* is the propagation time of the stress 
pulse through the aluminum cover. r is the 
propagation time of the stress pulse through the 
quartz guage. 



J. AZIMUTHAL MAGNETIC FIELD SYMMETRY 

To verify the azimuthal symmetry of the self-induced magnetic field 

O 

a series of pictures were taken at 9 = 5 intervals at (1.0, 6, 1.5). 
The magnetic probe was inserted through a holder in port five of the 
chamber (See Figure 4). The probe holder is circular and nests in a 
circular flange allowing rotation of the probe without disturbing the 
r or z coordinates. 

Figure 6l reproduces the probe response at 15 degree intervals from 
0 = n to 0 = tf/2. The gross aspects of the magnetic probe response 

are symmetric. 
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Figure 6l. rfegnetic Field Symmetry for (1.0, 8, 1.5). Scales are 
100 nsec and 100 mV per division. 



Unfortunately this experiment was conducted prior to discovery of 
the plasma expansion asymmetries. Therefore the oscilloscope scales 
were 100 nsec and 100 mV per cm. These scales and the z coordinate of 
1.5 cm do not allow resolution of the early magnetic field asymmetries 
expected with the plasma density anisotropy. Therefore, this 
experiment should be repeated with r = 0.5 cm, 0 = zero to n and 
z = 0.5 cm. The oscilloscope scales should be set at 100 mV and 
20 nsec per division. 



88 



IV. ANALYSIS OF EXPERIMENTAL RESULTS 



In this section two facets of the experimental results are examined. 

In section (A), a comparison is made between the particle content in 
the early plasma pulse and the number of atoms initially exposed to the 
laser radiation (to a depth equal to the skin depth). In section (B) 
the J X B force density is analyzed as a possible source for the radial 
confinement . 

A. ORIGIN OF THE EARLY PLASMA PULSE AND ITS PRESSURE DEPENDENCE 

Assuming a macroscopically neutral plasma and plasma symmetry about 
the axis of propagation, the total number of ions involved in the early 
plasma pulse is approximately 2 X lO 1 ^ ai ions. The skin depth for 
aluminum is 6l angstroms. The laser focal spot size is 0.15 cm. The 

O n 

total volume initially exposed to the laser radiation is 1.08 X 10“ cm. 
This corresponds to a total of 6.5 X 10^ aluminum atoms. This agrees 
to within one order of magnitude with the number of observed ions. If 
the early plasma pulse also has neutrals by the time it reaches the 
probe grid then this agreement will be even better. Therefore, at least 
in terms of total numbers, the early blowoff of plasma could be due to 
those atoms which are directly stimulated by the initial laser radiation. 
The pressure dependence of the early plasma has been previously 
described by Arifov [Ref. 5]* Assuming that the early blow off of 
material is independent of background gas density implies that there 
are initially 2.4 X 10^ ions/cm^ in the fast stream. Background pressures 

_5 

reach this density at 7 X 10 torr. This pressure is consistent with 
onset of the pressure dependent region. The fast streaming early plasma 
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further ionizes this background. However, when the background pressure 
is low, the additional ions do not add appreciably to the total. As 
background densities increase above the threshold there are 
significantly more ions created by the fast streaming plasma. As the 
background pressure is increased above 250 mtorr the fast streaming 
early plasma must be stopped in a short distance causing a reduction 
in early signal strength. 

B. ANALYSIS OF THE J XB FORCE DENSITY 

McKee [Ref. 38 ] using typical values for plasma expansion Into a 
250 mtorr background found that the J X B force density was much less 
than the radial pressure gradient force density. McKee also plotted 
this force density at 120 nsec for 250 mtorr of nitrogen background gas. 

He found that the magnitude of the largest force at this time is 

h O 

1.4 X 10 4 newtons/m. A reproduction of his computer produced plot is 
included as Figure 62. 

Bird [Ref. 8] has found magnetic fields of the order of 200 gauss 

at (0.4, TT/2, 0.75) for aluminum plasmas. This measurement was made at 

200 nsec after arrival of the laser pulse. The pressure was 0.1 microns 

which is close to the pressure independent region for both the magnetic 

field magnitude and magnitude of the early disturbance as reported by 

Arifov [Ref. 5]* The current corresponding to this 200 gauss field is 

2 

of the order of 400 amps/cm . Using maximum current and magnetic field 
densities indicates the presence of a J X B force density of approximately 
7 X 10^ nt/m^. This is consistent with McKee’s maximum values. 

At 500 nsec, the extent of the plasma (2 volt contour) is 
approximately r = 1.5 cm by z = 4.5 cm. Assuming syimetry with respect 
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Figure 62. J X B force density at 120 nsec for 250 mtorr of nitrogen 
background gas. The magnitude of the largest force at 
this time is 1.4 X 1(T nt/m3. 



to the target normal this corresponds to a total volume of 32 cm-3. If 
10 1 ? atoms are ejected and ten percent are ionized, then the ion density 

ill O 

would be on the order of 3 X 10 ions/cm. 

In section III. F. ion densities for the streaming main plasma of 
1.4 X 10 1 3 V p (V p is the probe response voltage) were found. Typical 
probe response voltages are on the order of ten volts for the main plasma. 
This corresponds to a total ion density of 1.4 X 10^ ions per cm^. 

These two results are in good agreement. The total plasma density 
including neutrals is then taken to be 10^ cm - 3. 
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Assuming a plasma temperature of three electron volts, and a 
characteristic density gradient length of / = 0.8 cm, the pressure 
gradient force density becomes: 

(17) nkT/ f = 10^ newtons/m^ 

Using either the maximum J X B force density found by McKee or that 
found on page 90 we see that the J X B force density is of the same 
order of magnitude as the pressure force density found in equation (17)* 
Therefore the J X B force density is capable of at least partial 
confinement of the plasma and in any case should be taken into account 
in describing the dynamics of the expanding laser produced plasma. 

In order for this confinement to continue for the observed times, 
the characteristic diffusion times must be at least of the same order 
as the observed effect. Otherwise the plasma will diffuse through the 
magnetic field. 

The magnetic field diffusion time is given by: 

( 18 ) r a = K ot 

The characteristic length ( / ) of the magnetic field at 200 nsec is 
estimated from McKee's data to be 0.8 cm. The D.C. conductivity of the 
plasma is given by: 

(19) ° = n e e 2 /nfe v e 

where v e is the electron collision frequency. For electron-ion and 
electron-neutron collisions the collision frequencies are given by 
Tanenbaum [Ref. 46]. 
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(20) w el 



= 3-62 X 10“ 6 n;L T e 3/2 InA sec 1 

(21) v = 2.60 X lO^n T 1/2 a 2 sec _1 
en n e 

Where T is in °K, n in particles/m 3 , In a is typically about ten, a is 
of the order 10“ 10 m. Assuming T e = 4 eV, nj_ = 10 20 m -3 and = 10 21 m“ 3 
the collision frequencies are found to be = 3.62 X 10 3 sec -1 and 
u en = 5*6 X 10 7 sec Substituting these values into equation (19) 
reveals a D.C. conductivity of 7-8 X 10 3 mho/m. Using this value in 
equation (18) we find that the magnetic field diffusion time is of the 
order: 



(22) r m = 323 nsec 

This value appears consistent with McKee's work and explains how the 
plasma could be radially confined within the time scale surveyed in this 
investigation (t< 800 nsec). While not plotted, an examination of the 
data indicates that the plasma becomes very diffuse after 900 nsec. 

Unis is due in part to diffusion but is more probably due to general 
cooling and recombination. 



93 



V. SUMMARY AND CONCLUSIONS 



Ihe early disturbances have been separated into "reverse" photo- 
electric response due to photoemission of electrons from the probe and 
a series of at least two additional disturbances traveling with constant 
speeds of 1.1 X 10® cm/sec and 5-9 X 10? cm/sec. The first of these 
disturbances has been analyzed in detail. 

Ihe early disturbance is found to propagate along a line describing 
the path of the reflected laser pulse. The disturbance is symmetric 
about this propagation axis. 

After separating the photoelectric response, consideration must be 
given to the possibility that the remaining disturbance is created by an 
ionizing potential wave or a blast wave. References 2, 9, 22 and 47 
provide some insight into the ionizing potential wave phenomenon. In 
this experiment, the background gas density is approximately 9 X 10-*--*- cm - ®. 
Assuming that the temperature of the background is three electron volts 
after ionization by the potential wave we can use equation (25) from 
Appendix A to obtain an estimate of the quiescent plasma ion density 
necessary to generate the observed probe current of 0.66 amperes at 
( 0 , 0 , 1 . 0 ). 

(26) n. = 1.67 I + ( rrij/kT e ) 1/2 / eAp (cm" 3 ) 

Ihe average molecular weight of air is m = 28.96 AMU = mj_. The surface 
area of the electrodes is 0.0146 cnf*. Inserting these values into 
equation (25) indicates that the required ion density is 1.5 X 10-*-^ ion 
charge s/cm®. This ion density is three orders of magnitude greater than 
the background gas density. Ibis argument appears sufficient to exclude 
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ionization of the background gas by any means as a possible source for 
the early disturbance. 

The higher densities observed can only be supplied from the target 
in the form of a fast blowoff of material. This possibility has been 
considered by many researchers [Refs. 10, 15, and 24]. However, a 
review of the literature indicates that this early blowoff is generally 
thought to be composed entirely of electrons. The results of the present 
research in which transverse electric and magnetic fields were introduced 
indicate that the early disturbance is composed of ions and electrons. 

Ion densities of 2.4 X 10^ cm“3 have been measured in this fast streaming 
plasma. Electron densities may be slightly higher than this creating a 
macroscopically negative plasma potential. This would explain the 
previous observations. 

The pressure dependence of the early plasma has been previously 
described by Arifov [Ref. 5]. Assuming that the early blowoff of material 
is independent of background gas density implies that there are initially 
2.4 X lO-*-^ ions/cm^ in the fast stream. Background pressures reach this 
density at 7 X 10 ^ torr. This pressure is consistent with onset of the 
pressure dependent region. The fast streaming early plasma further 
ionizes the background. However, when the background pressure is low, 
the additional ions do not add appreciably to the total. As background 
densities increase above the threshold there are significantly more ions 
created by the fast streaming plasma. As the background pressure is 
increased above 250 mtorr, the fast streaming early plasma must be 
stopped in a short distance causing a reduction in early signal strength. 
In this investigation the early plasma density is strongly attenuated 
between 0.8 and 2.0 centimeters. However, a small sigial is still 
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observable out to the limits of this survey (5 cm). Ftesidual 
ionization decreases rapidly so that by 100 nsec there is little 
remaining evidence of the passage of this early pulse. 

Ihe main plasma speed has been measured at 1.1 X 10? cm/sec. Early 
time anisotropies and asymretries are observed in the expansion of the 
main plasma. This results in an initial high density area on the side 
of the target away from the incident laser pulse. Ihis asymmetry is 
short lived and for times greater than 120 nsec, the plasma can be 
described as symmetric about the target normal. This late time 
symmetry in the free expansion is evident even when the plasma is created 
in a deep hole and gains momentum in a direction off the target normal 
before reaching the surface of the hole. In this investigation a one 
centimeter deep hole was excavated into an aluminum alloy target. 
Throughout the excavation there was little observed change in the plasma 
expansion. This suggests that the expansion dynamics are independent 
of the manner in which the plasma reaches the surface of the target . 

Previous investigators [Ref. 23] have observed approximately 10-^ 
particles ejected per incident laser pulse under similar conditions. 

In this hole drilling exercise approximately lO"^ atoms were ejected per 
incident laser pulse. Most of this material was neutral and completely 
coated the exposed probes. The reason for identifying the additional 
material as neutral is that little change was observed in the ion 
densities at probe one and probe two though heavily coated with aluminumj 
never did register a significant ion density. 

Examination of the cross sectioned crater indicates that a 
significant amount of the target damage is generated by dynamic shock 
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waves rather than by thermal processes. This is especially true of the 
base of the crater where whole chunks of material are being ejected, 
possibly by dynamic shock processes. 

Hie main laser plasma is observed to expand very rapidly in the 

radial direction to a diameter of approximately 2.4 cm. This maximum 

radi a l extent is achieved at an axial distance of 0.3 cm. At this point 

the plasma plume ceases to grow in the radial direction and expands 

symmetrically along the target normal with a cylindrical shape. Peak 

plasma densities in the main plasma pulse are approximately 
14 

2.5 X 10 cm . Radial density gradients remain high for at least one 
microsecond while the axial density gradient decreases steadily with 
increasing time. This behavior is strongly suggestive of the presence 
of a radial confining force acting on the expanding plasma. The self- 
induced magnetic fields generate a J X B force density which is seen to 
be comparable (within an order of magnitude) with the radial pressure 
gradient force density. This confining force is postulated as an 
important element in creating the observed cylindrical shape of the 
laser produced plasma. This J X B force must be included in any dynamic 
model of the expanding plasma. 

Table (3) summarizes the properties of the various plasma features. 



PLASMA PEATUKE 


SPEED 
(cm/ sec) 


ION CHARGE 
DENSITY (cm-3) 


TOTAL 

NUMBERS 


TOTAL' ENERGY 
(joules) 


Photons 


3.0 X 10 10 





_ 


8 X 10 -5 


First Early Plasma 


1.1 X 10 8 


2.4 X 10 12 


2 X 10 1 3 


0.55 


Second Early Plasma 


5.9 x 10 7 


3.0 X 10 12 


1 X 10 13 


0.079 


Main Plasma 


1.1 X 10 7 


2.5 X lO 1 ^ 


1 X 10 16 


1.1 


Target Thermal Energy 


— 


— 


- 


1.5 
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Hie total energy radiated by photons with wavelengths between 
0.12/x and 0.267 m has been previously computed. The energy content in 
each of the two early plasmas assumes that each of the plasma particles 
has a directed speed equal to the speed of the plasma front. The 
ionization energy and assumed early plasma temperature of ten electron 
volts contribute a negligible amount of energy. 

The energy of the main plasma is much more difficult to compute. 
Using the methods outlined by Dawson [Ref. 14] the temperature in the 
main plasma body is estimated to be T q = Z + 1). Early in the 

expansion, the average ionic charge is assumed to be Z = 3- This gives 
a plasma temperature of 170 eV and a total plasma energy of 1.1 joules. 

Examination of Figures 37 to 44 indicates that only a small portion 
of the particles in the main plasma have the directed velocity seen in 
the plasma front. The bulk of the main plasma does not move from the 
immediate vicinity of the target . The sum of those energies accounted 
for in Table (3) is 3*2 joules or about 65 percent of the energy in the 
incident laser beam. The energy absorbed in the target by conduction 
and shock wave absorption has been estimated at 20 to 30 percent of the 
incident laser pulse energy by Peter Krehl at the Stanford Research 
Institute (personal comrmunication) . 

This account of the incident energy neglects energy contained in the 
neutral particles which constitute 90 percent of the ejected material, 
reflected laser pulse energy and plasma radiation at frequencies other 
than those between 0.12 and 0.267 microns. These mechanisms must 
account for the additional 35 percent of the incident energy . 
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VI. AREAS WARRANTING FURTHER RESEARCH 



Because of the early time asymmetry of the expanding laser plasma, 
a three-dimensional map of both the early and main pulses is necessary 
to fully understand their extent. In conjunction with this, the data 
obtained in this investigation should be used to provide a computer 
generated map of the plasma density. 

Development of plasma temperature diagnostic techniques at NPS would 

allow production of similar maps for the plasma's isotherms. With this 

information, the school’s computer should be used to plot the cross 

product of density and temperature gradients. This would predict the 

magnitude and direction of the magnetic field due to the VkT e X Vn e 

e 

source term described by McKee and Bird. These results could then be 
compared with McKee's measurements [Ref. 38 ] in the pressure independent 
region described by Bird [Ref. 8]. 

This investigation has revealed asyirmetries in the early expansion 
of the laser produced plasma from a plane target. The dependence of 
these asyirmetries on the angle of incidence of the laser pulse requires 
investigation. In conjunction with this, an accurate measure of the 
reflected laser pulse energy at 1.06 microns should be correlated with 
the angle of incidence. 

The pressure dependence of these early time asymmetries could help 
explain the pressure dependence of the early plasma density and 
magiitude of the self-induced magnetic field strength [Refs. 5 and 8]. 
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Measurement of currents in the target generated by the inhomogeneous 
blowoff of ions and electrons and by the self-induced magnetic fields 
could help resolve the induced magnetic field and plasma density 
characteristics at distances too close to the solid target material to 
be measured with external probes. 
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APPENDIX A 



I. DOUBLE PROBE THEORY 

When a potential difference is applied between two electrodes 
immersed in an ionized but macros copically neutral plasma, the positive 
ions are attracted to the negative electrode and the plasma's electrons 
are attracted to the positive electrode. This action tends to destroy 
the plasma's neutrality by creating space charge electric fields. The 
plasma responds -to this perturbation by creating a sheath which restricts 
the charge separation to a short distance from the intruding electrodes. 
Ihis sheath absorbs the potential difference in a very short distance 
which is of the order of a few Debye lengths. The Debye length is equal 
to: 

(23) A =69 X (T/n ) 1/2 (meters) 

D e 

14 

Assuming a temperature of 100 eV and an electron density of 10 

-3 -4 

cm , the Debye length is 7.4 X 10 cm. The bulk of the 

plasma remains charge neutral. 

When the surface area of the two electrodes is equal and the entire 
system (electrodes, biasing power supply, scopes, etc.) is isolated from 
ground it is said to "float" with the plasma potential. 

If the biasing potential is initially zero and the double probe is 
iimiersed in a neutral plasma with T. = T then the random electron flux 
to the probe tips will be larger than the random ion flux due to the 
higher electronic thermal velocity. This will cause the probe tips to 
collect a net negative charge until the potential change is just 
sufficient to allow an equal flux of ions and electrons. The probe tips 
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will then be at a potential below that of the plasma. This potential is 
called the floating potential. This is illustrated in Figure 63- 




Electrode 1 Electrode 2 

Figure 63 . Probe potential with zero bias applied. 

If the number two probe is now biased slightly negative, V 2 will 
become more negative and will become slightly less negative. Thus 
more electrons will flow to probe one. This results in a positive 
current flow to probe two. When the number two probe potential becomes 
negative enough to exclude electrons and collect the saturation ion 
current, probe number one will be at a potential just sufficient to 
allow collection of the net electron flux necessary to cancel the ion 
current to probe one. This condiction satisfies Kirchoff's Law: 
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The potentials involved in this situation are depicted in Figure 64 and 
the current flows in Figure 65 . 

An important advantage of the double probe is that the current in 
the system is limited by the ion saturation current. In the absence or 
imbalance of ions, the potential of the entire system changes so that the 
maximum current drawn is always limited by the ion densities present. 
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Figure 64. Collection of saturation ion current at probe two. 




Figure 65 . Currents flowing in a double probe collecting saturation 
ion current at probe number two. 



For a macros copically static plasma, the saturation ion flux to a 
cylindrical probe of surface area (A) is given by [Ref. 45]: 

(25) n ± = 1.67I + (m i AT e ) 1 / 2 /eA p (cm"3) 
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The electron temperature (T ) can either be estimated or found from 
the slope of the probe characteristic at the origin [Ref. 11). 

(26) dl/dV =_e_.(i 1+ ) . (i ) 



If the areas of the two electrodes are equal and the product of the 

probe bias voltage and the average ionic charge is much less than the 

average thermal energy of the ions, then i 1+ = i 2+ . Substituting this 

into equation (26) and solving for T we find that™ 

e 

(27) T = el. (l/dI/dV)/2k 

It should be remembered that these equations are valid only in a 
static plasma. In the current investigation, the initial random thermal 
energy is rapidly converted to ordered flow kinetic energy. This is at 
least true for plasma expansion into near vacuum conditions. Dawson 
[Ref. 14] finds that the expansion appears to have 'taken place by the 
time the plasma has expanded to dimensions of 0.01. to 0.1 cm. If the 
flowing plasma is to remain charge neutral, then n e = n^Zj_ and the flow 
velocities of ions and electrons must be equal. Therefore , the ordered 
kinetic energy of the plasma is predominantly due tso the ions because 
of their greater mass. 

Assuming that T 0 = Tj_ = 100 eV and a flow velocity of 10? cm/sec, 
the ratio of ion thermal to ion directed speed is 2.5 X 10“^ and the 
ratio of the electron thermal speed to flow speed is approximately 13. 
Therefore the smaller current will be that of the Lons. 

Now I + = V/R where R is the resistance across which the current is 
measured. This is also the average number of positive charges collected 
by the probe per unit time. This can be expressed as: 
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(28) V/R - n^AVy, 

Where rij_Z e is the average electronic charge per unit volume and (A) is 
the cross sectional area of the probe perpendicular to the plasma stream 
velocity (V^). 

Equation (28) does not include secondary emission at the probe tips 
due to ion bombardment. Koopman [Ref. 30] has measured this effect and 
finds that the secondary emission coefficient varies between zero and 
0.5 electrons emitted per incident ion. He also found that secondary 
emission was not strongly dependent on the ion kinetic energy in the 
range of energies encountered here. Therefore, the secondary emission 
can be treated as a constant. To include this effect, equation (28) 
can be rewritten to include an average number of secondary emissions 
per incident ion: / 




(29) V/R = n ± Z e (l + O A V 



Since A, V, R and V^p are measured quantities and ( is a constant, 
equation (29) provides a reasonable estimate of the absolute plasma 
density. For equation (29) to be valid we require that the probe 
electrodes be of equal area and that the plasma be homogeneous across 
the probe tips. The probes were carefully constructed to satisfy the 
first condition. Previous workers at NPS have been careful to keep 
the plane containing the probe tips perpendicular to the direction of 
the streaming plasma. However, the probe tips are about one millimeter 
apart and the characteristic density gradients have lengths of five 
millimeters. Therefore, it was suspected that this orientation was not 
critical. To test this hypothesis, the probe was maintained in the 
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same geographical position and rotated through a variety of angles. 
Figure 60 shows the results as recorded from the oscilloscopic tracings. 
Ihese reproductions show that the signals are essentially unaffected by 
rotation of the plane containing the probe tips. 

Qi the basis of this experiment, the plane containing the probe 
tips was maintained perpendicular to the z-axis throughout the 
experiments regardless of the expected plasma flow velocity. 




Figure 66. Effect of rotating the plane containing the probe tips with 
respect to the targets z-axis. The angles (0) given are 
the angles between the plane defining the probe tips and 
the target normal, (a) recorded at two volts and 100 nsec 
per division and (b) is recorded at 500 mV and 20 nsec per 
division . 
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II. DOUBLE PROBE CONSTRUCTION 



Figure 67 depicts the probe circuit used throughout this investigation. 
In one instance the four microfarad capacitor was replaced with a 0.5 
microfarad capacitor. No detectable difference in the probe response was 
noted. 




Figure 67. Double Probe Circuit 



The following tests have proven effective to Insure proper probe 
operation. 

!• Check each electrode for continuity (no resistance between 
(A) and (I) and between (B) and (J)). 

2. Insure isolation of the probe circuit from the aluminum housing. 

( >100 kfl between (H) and (E) and between (C) and (E)). 

3. Insure adequate grounding of the stainless steel probe shield, 
(zero resistance (F) to (E)). 

4. Check the current limiting resistor. (>5kQ(G) to (B)) 

5. Check the current viewing resistor for one ohm. ((D) to (C)). 

6. Check -the circuit for continuity. (No resistance (H) to (D) 
and one ohm from (H) to (C)). 

lYbst of the glass enclosed probes were made by a professional glass 
blower. These probes had very fine tungsten electrodes with a diameter 
of 0.013 cm. The electrodes were 0.25 cm long and were set approximately 
0.08 cm apart. The diameter of the glass tip was approximately 0.1 cm. 

At one point in the investigation, a larger signal to noise ratio 
was desired. This was necessary when examining the early disturbance at 
distances greater than two centimeters from the target. In order to 
achieve this higher ratio a double probe was constructed with large 
electrodes. Production of this probe required little technical expertise 
and the rret hod is repeated here as it requires a minimum of equipment and 
could be duplicated in any laboratory. 

The electrodes were constructed from 0.05 cm diameter tungsten rods. 
These rods were cut to approximately four centimeter lengths and one end 
of each rod was spot welded to a 0.05 by 0.3 cm nickel interface. A 
solid copper lead was spot welded to the other end of the nickel interface. 
This is depicted in Figure 68. 
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Copper Wire 



Nickel 

Interface’ 



Figure 68. Tungsten electrode connection to the copper lead through 
a nickel interface. 



Two of these electrodes were then cut to lengths differing by one cm 
(2.5 cm and 3*5 cm). It should be noted that tungsten has a tendency to 
split when snipped. The most successful cutting method involves notching 
the tungsten rod with a file and then breaking it over a sharp comer 
much as glass rod or tubing is cut. The tips were then etched in a 
saturated solution of sodium hydroxide at the same time. The etching 
process is depicted in Figure 69 . 





Carbon 

Electrode 



O- 



~|+ AC Power 
Supply 

J- 



Sodium Hydroxide 

Saturated Solution 



.Tungsten 

Probes 



Figure 69 . Etching the tungsten electrode tips. 



The etching time was ten minutes and produced identical tips as depicted 
in Figure 70. 
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ip diameter .014 cm 
base diameter .034 cm 



protruding length .1905 cm x 

Figure 70. Etched Tungsten Electrode 

In the next step, a five millimeter diameter pyrex tube was cut to 
a length of 20 cm. Two 0.75 mm outside diameter, 0.60 mm inside 
diameter alumina tubes were cut to lengths of 1.5 and 2.5 cm and a piece 
of RG-17VU microcoaxial cable was cut to a length of 30 cm. One end 
of the coax was stripped and inserted inside the stainless steel shield 
which is then inserted inside the pyrex tube. At this point the 
electrodes are soldered to the microaxial cable. Figure 71 depicts the 
probe in this state of completion. 







To insure proper insulation, all exposed leads in the two electrode 
circuits are covered with a thin coat of "RTV. " When this is dry the 
probe should be checked for circuit continuity, (zero ohms (B) to B') 
and between (A) and (A’)* 
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Ihe supporting structure for the electrode tips is made of liquid 
procelain. The exposed electrode tips are lightly tied together with 
thread and the entire circuit is stretched between two clamps . A thin 
coat of procelain is then applied. After the first coat is dry, the 
probe may be removed from the clamps. Several additional thin coats 
should be applied until the desired strength is achieved. Figure 72 
depicts the finished probe. 



LIQUID PORCELAIN COVERING 








Figure 72. Finished Double Probe (a) internal configuration (b) 

external dimensions . . * Insure that the liquid procelain 

is drawn into the aparatus here to provide adequate 
strength. 



'Because of the larger size of this probe, there was some concern 
that the plasma might be perturbed sufficiently to distort the data. 

To determine what effect this larger probe had on the plasma features a 
comparison was made between the probe response of the finest probe 
available and this slightly larger probe. The comparison is duplicated 
as Figure 73* As can be seen, all features are clearly retained, the 
primary difference being a larger signal from the probe with the larger 
electrodes. Because of the additional size of this probe it was not 
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used in any experiments where additional diagnostic information was 
taken downstream. 




Figure 73- Comparison of the density profile sampled at (0, 0, 1.0) by 
(a) the large probe biased at -20 volts and (b) by a very 
fine probe biased at -15 volts. 
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APPENDIX B 



I. PRESSURE PROBE THEORY 

When a force is exerted on a piezoelectric crystal such as quartz, 
it results in the formation of a stress within the 'crystal. This stress 
distorts the ionic structure of the crystal resulting in an increase in 
the polarization of the sample. This process is depicted in Figure 7^. 

STRESS 
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♦0.0.0 

P+AP _ „ 


o o o 




"o o o 




1 t 



STRESS 



(a) (b) 

Figure (a) Unstressed Piezoelectric Crystal 

(b) Stressed Piezoelectric Crystal 

The stress changes the polarization by AP, the induced 
piezoelectric polarization. 



Recall that the electric displacement vector is given by: 

(30) "D = (E +T~ 
x xx 

For a one dimensional strain, the piezoelectric equation is: 



(31) D = Td + (E 
x 21 

When D x = electric displacement vector, T = stress, d 2 ^ is the 
Piezoelectric Strain Constant, E = electric field and c is the 
permittivity . 
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For an X-cut quartz crystal [Ref. 18] 

( a ) a = 2.25 X 10" 12 coulombs/newton 

(b) ( = 4.06 X 10 -11 farad/meter 

In the static case the open-circuit voltage is given by [Ref. 4]. 

(32) g = d /c = 0.055 fF/tw volt-meters/newton 
and the short-circuited charge for a given applied force is: 

(33) Q = d 21 Ff/t 

where F is the force which is considered positive for an extensional 
stress, i = crystal length, w = crystal width and t = the crystal 

thickness. 

Die material presented above is satisfactory for low frequency impulses. 
In this investigation, the anlysis is fisher complicated because the 
r i„« produces a stress wave with a short wave length. Graham, et al., 
[Ref. 18] have analyzed the short-circuit piezoelectric current induced 
by the application of a rapidly changing impulsive load. The load 
propagates as a stress wave along the x-axis of a quartz disk. The 
displacement current generated by this stress wave will be: 

(34) i = A dD/dt 

D 

Where (A) is the electrode! area and D is the displacement vector. 

Substituting equation (3D for (D) we find that. 

i i f 

(35) /D(x)dx = /p(x) dx + e/E(x) dx 
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Where l is the thickness of the quartz disk. By assuming that a short 



circuit exists between the electrodes and that the crystal has constant 
permitivity the last integral in equation (35) becomes zero [Ref. 18]. 



The conductivity of the crystal is essentially zero so that equation 

(35) becomes: 



At this point it is assumed that the piezoelectric polarization, 
P(x), is proportional to the x-component of the stress by a coefficient 
(f). This coefficient is independent of time and stress for a given 
stress range. Therefore: 

(38) P(X) = f o (x) 

Where a (x) is the x-component of the stress. This assumption is 
verified by Graham [Ref. 19]. With the assumption that the quartz is 
linearly elastic we can write: 

(39) a (x,t) = o (x - U s t) 

Where U_ is the wave propagation speed in the crystal. 

S 

The solution for the current from equations (3*0 and (37) is there- 
fore found to be: 



(36) Jt E(x)dx = 0 



0 



(37) 1 f P(x)dx = D 

! 0 



(40) i = AdD/dt 
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Wh ere o 0 is the x-corrponent of stress at the stress input electrode 
and is the x-conponent of stress at the rear electrode. 

For early times in the wave transit, o =0 and equation (40) 
becomes : 



(HD i = Ms", 

l 

Equation (4l) predicts that for times less than the wave transit time, 
the current is proportional to the stress which is proportional to the 
force applied to the crystal which in the present investigation equals 
the momentum transported to the crystal by the plasma. Therefore in 
this investigation the pressure probe response is directly proportional 
to the momentum transport in the plasma. 
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II. PRESSURE PROBE CONSTRUCTION 



Stanford Research Institute provided the pressure probe used in 
this investigation. It is a Sandia Quartz Gauge. The time resolution 
of this type gauge is on the order of a few nanoseconds [Ref. 18]. 
Figure 75 depicts the construction of the gauge. 




Figure 75* Sandia Quartz Gauge. 

This gauge used a guard ring configuration described in detail by 
Graham [Ref. 18]. The configuration is pictured in Figure 76. 




Figure 76. Guard Ring Electrode Configuration 

The inner portion of the quartz disc is isolated electrically from the 
outer portion by separating the vapor coated electrode into two regions 
as shown above. 
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This configuration helps to eliminate field distortions. At the 
outer edge of the disk, the discontinuity in electric potential and 
dielectric permittivity causes electric field fringing similar to that 
found in parallel plate capacitors. The guard ring eliminates this 
distortion by restricting observation to the central region of the disk 
where there is negligible fringing. 

The guard ring also helps to eliminate distortions due to the 
generation of unloading waves. As the stress pulse moves through the 
disk in the axial direction, boundary conditions generate shear and 
dilatational waves immediately behind the wavefront at the lateral edge 
of the disk. These unloading waves then propagate laterally inward 
from the lateral boundary and the crystal is no longer in a state of 
one dimensional strain. The central region of the disk will be in one 
dimensional strain only for the first wave transit time if the width of 
the outer electrode is such that an unloading wave does not reach the 
central region during the first wave transit time. 
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